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1 Introduction 

The purpose of this paper is to give a complete derivation of the hmiting distribution of large 
Frobenius numbers outlined in [T] and fill some gaps formulated there as hypotheses. We start with 
the basic definitions and descriptions of some results. 

Consider n mutually coprime positive integers ai, a2, ■ ■ ■ , an. This means that there is no r > 1 
such that each aj, 1 ^ j ^ n, is divisible by r. Take which later will tend to infinity and will be our 
main large parameter. Introduce the ensemble Qn of mutually coprime a = (ai, . . . , an), 1 ^ aj ^ N, 
^ ^ j ^ n and Pjy be the uniform probability distribution on Qn. For each a £ Qn denote by F{a) 
the largest integer number that is not representable in the form x = xiai + ■ ■ ■ + Xnan, where Xj 
are non- negative integers. F{a) can be considered as a random variable defined on Qtv- The basic 
problem which will be discussed in this paper is the existence and the form of the limiting distribution 

for the normalized Frobenius numbers /(a) = —F{a). The reason for this normalization will 

be explained below. 

The case of n = 2 is simple in view of the classical result of Sylvester (see [7] ) according to which 
F{ai, 02) = aia2 — ai — 02. It shows that in a typical situation F grow as N'^. The first non-trivial 
case is n = 3 where F{a) grow as N^/"^ It is known (see [10]) that the numbers F{ai,a2,a3) have 
weak asymptotics: 



For arbitrary n the only result known to us is the following theorem proven in [T]. 
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Theorem 1. Under some additional technical condition (see 11]) the family of probability dis- 
tributions of fi^{a) = ^"T— -F(a) is weakly compact. This means that for every e > one can find 

V = V{e) such that 

I J 1 F{a) ^ P j ^ 1 - e. 

In this theorem e,T) do not depend on N. It also imphes the existence of the hmiting points 
(in the sense of weak convergence) for the sequence of probabihty distributions of /7v(a)- As was 
aheady mentioned, in this paper we shall study the limiting distribution of /Ar(a) = j^j2F{a), 
a = (ai, 02, 03) as ^ cxD. This distribution is not universal and will be described below. 

Take any p, < p < 1, and consider its expansion into continued fraction 

p=[hi,h2,...,hs,...] (1) 
where hj ^ 1 are integers. If p is rational then the continued fraction ([1]) is finite. The finite 

Ps 

continued fractions p = [hi, . . . ,hs] = — are called the s-approximants of p. The numbers Qs satisfy 
the recurrent relations 

Qs = hsQs-l + qs~2, s^2 (2) 

Introduce the Gauss measure on [0, 1] given by the density 7r(x) = i^2{i+x) • "^^^^ the elements of the 
continued fraction ([1]) become random variables. It is well-known that their probability distributions 
are stationary in the sense that the distributions of any hm-k, ^m-fc+i) • • • > ^m- • • • 1 ^m+fe do not 
depend on m. We shall need the values of s = si, such that q^i is the first qg greater than \/iV. It 
was proven in [5] that qs^/\fN have a limiting distribution. More precisely, the following theorem 
is true. 



Theorem 2. Let k be fixed and s{R) be the smallest s for which qg ^ R. As R ^ 00 there exists 
the joint limiting probability distribution of hg(^ji'^_}., . . . , 

In the paper |llj the analytic form of this distribution was given. 

Consider the subensemble Q'f^ C Qn for which 01,03 are coprime. Then there exists a^^(mod 

03), 1 ^ a^^ < 03. Denote p = The expansion of p into continued fraction will be needed 

below. Clearly, p is a rational number. However, the following theorem is valid. 

Theorem 3. As before, consider si such that qsi-i < \/N < g^-^. Then in the ensemble Q^"* 
equipped with the uniform measure, for any k > and N ^ 00 there exists the joint limiting 
probability distributions of hs-^-k, ■ ■ ■ , hsi+k which coincides with the distribution in theorem 2. 

A stronger version of theorem [3] is also valid. 
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Theorem 4. Let the first elements of the continued fraction for p be fixed: /ii, /i2, • • • ,hi. Then 
under this condition and as N ^ oo the conditional distributions of ^j=, hg-i^-k, ■ ■ ■ , hg^^-^k converge 
to the same limit as in theorems and O 

All these theorems will be proven in section [3l Now we can formulate the main result of this 
paper. 

Theorem 5. There exists the limiting distribution of f]\f{a) = fN {{0,1,02,03)), (01,02,03) € Qn 
as N ^ 00. 

The proof of the main theorem is given in section [2j First we consider the ensemble Q^^^ and 
then explain how to extend the proof to Q^. 
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2 The limiting Distribution of /iv(a). 



Return back to the case of arbitrary n. Introduce arithmetic progressions 



Hr = {r + man, m ^ 0}, ^ r < O; 



For non-negative integers xi, . . . , Xn-i such that xioi + X2O2 + • • • + 2;„_io„_i G we write 



xioi H h x„_io„_i = r + m(xi, . . . , x„_i)on. 



Define m(r) = min m(xi, . . . , x„_i) and put 

X\...,Xn-l 



Fi{a) = max min (r + m(xi, . . . , Xn-i)on) = 

xiai^ hx„_ia,^_ieYlr 

= max min {xiOi + . . . + a„_io„_i) . 

0^r<a„ xiai+...+Xn+a„-i = r{mod a„) 

It was proven in [3] that F{a) = Fi{a) — o^- A slightly weaker statement can be found in [1]. Since 
in a typical situation aj grow as N while -Fi(o) grow as N^~^"^ (see also |T]) the limiting behavior 
of — ^^ ^\ and '^^ ^"1 is the same, but the analysis of ^'^^^l is slightly simpler. Let us write for 
n = 3 

xioi + X2O2 = r + m(xi, X2)a3 



or 

xiOi + X2O2 = r(mod 03) (3) 
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We assume that 01,03 and 02,03 are coprime. Therefore there exists o^^ ^, 1 ^ o^ ^ < 03, such that 
oi • o^^ = l(mod 03). Choose a^^ so that 1 ^ a^^ < 03 and rewrite ([3]) as follows 

xi + 012x2 = ri (mod 03) (4) 

where 012 = aj^^a2(mod 03), < 012 < 03 and ri = ro|f ^(mod 03), ^ ri < 03. From ^ 

012X2 = (ri - xi)(mod 03) (5) 

The expression ([5]) has a nice geometric interpretation. Consider S = [0, 1, . . . , 03 — 1] as a "discrete 
circle". Let TZ be the rotation of this circle by 012, i.e. 

TZx = X + oi2(moda3). Then TVx = x + pai2(modo3) and ([5]) means that ri — xi belongs to 
the orbit of under the action of TZ. From the definition of -Fi (o) 

Fi(o) = max min (xioi + X2O2) = 

0^r<a3 xiai+X2a2=r( mod 03) 

At3/2 ■ / Xl 01 X2 02\ , , 

= Jy ' max mm —= 1 — — (6) 

0sSri<a3xi+X2ai2=ri {mod aa) \^/N N ^/N N J 



Choose /i*^-'^ = {h!(\ . . . ,hm), j = 1,2,3 and denote by Q^^f^ii) f^(2) 1^(3) ensemble of o = 
(01,02,03) e such that the first m elements of the continued fractions of ^ are given by 

j = 1,2,3. This step means the localization of the ensemble Q^^^ ■ It is easy to see that for 
every e > one can find rational ai, 02, 03 and m such that |^ — aj| ^ e, 1 ^ j ^ 3. Then in ^ 
one can replace ^ by aj. Since -^L will take the values 0(1) the whole expression in ([6]) takes 
values 0(1) and instead of ([6]) we may consider 

f Xl X2 \ , . 

max mm I ai H — -1=^2 I (' j 

ri xi+ai2X2=ri (mod 03 ) \ V V / 

with the error 0{e). We assume that in the ensemble Q^^\(i) ^(2) ^(3) 'we also have the uniform 
distribution. 

We shall need some facts from the theory of rotations of the circle. According to our assumption 
oi2 and 03 are coprime. Therefore TZ is ergodic in the sense that 7?."^ = Id and 03 is the smallest 
number with this property. Put p = ^ and write down the expansion of p into continued fraction: 

p = [hi, /i2, • • • , hsg]. Let also be ps = [hi, /i2, ■ ■ ■ , hg] = ^ and si is such that Qs^-i < VN < qs^^. 

It will be more convenient to consider the usual unit circle instead of S and use the same letter 
TZ for the rotation of the unit circle by p. Introduce the interval A^J'^ bounded by and {qpp} and 
^(p) ^ 7^i/\(p). Using the induction one can show that Af\ ^ j < Qp+i and Af'^^\ ^ j' < qp 
are pair-wise disjoint and their union is the whole circle except the boundary points (see [5j). Denote 
by r/(p) the partition of the unit circle into A^j'\ A^J^~^^\ Then rj'^P^^') ^ r](p) in the sense that each 

element of rj^^^ consists of several elements of rj^P'^^'^ . More precisely, A^ consists of hp elements 
A'-f and one elements A^~^^\ The partitions r/^^^ show how the orbit of fills the circle. 
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Return back to the discrete circle S. The partitions rj^^ can be constructed in the same way as 
in the continuous case. We have to analyze 

max mm —=ai H :=ct2 (o) 



X'i+a\2X2=i'i{ mod az) 

for given ai, 02, < a\, < 1. 

Lemma 1. There exists some number Ci (01,02) = Ci such that for any ri the point xi giving 
min (^-^^ OL\ + a2^ under the condition is such that t\—xi fxi +012^2 = ri (mod 03) ) is an end- 
point of some element of the partition r]^^^~^'^^\ Here mi > is such that qsi + mi/qs^ < Ci{ai,a2) 

The proof is simple. In any case ri — xi is an end-point of some element of the partition 7^(^*1+™-!). 
If mi is too big then is too big because it takes too much time to reach an end-point of 
which is not an end-point of one of the previous partitions. We can choose yi so that ri — yi will 
be an end-point of some element of and the linear combination -^Cii + ;^^^2 is smaller. This 
completes the proof of the lemma. 

Its meaning is the following. If ri — xi is an end-point of ry(^i+'"i) with too big mi then X2 will 
be also too big. 

Lemma 2 shows that xi also cannot be too big. 

Lemma 2. There exists an integer m2 > depending on 01,02 the ratio qsi/N and the elements 
of the continued fraction hg^, hg^+i, . . . , hs^+m2 of p such that for any ri the interval [ri — xi, ri] 
corresponding to the minimum of 

Xl X2 

Oil H 1= OL2 



has not more than m2 elements of rj^^^K 

The proof is also simple. If the number in question is too big then will be too big. Therefore 
for given ri min can be attained at a point which is closer to ri. 



The values of qs^/vN and hs^,hs^+i . ■ ■■, hsi+m2 determine the structure of the partitions 

T)(*l) rt{^l+m2) 

'/ 5 ■ ■ ■ ) '/ 

The conclusion which follows from both lemmas is that for each ri we check only finitely many 
Xl and X2 and find min(xiOi -|- a;202) among them. The number of points which have to be checked 
depends on 01,02, ^ and hg^,.. .,hg^+rn2- 

Now we remark that ri must be also an end-point of some element of the partition Indeed, 

if ri increases within some element of rj^^^^ then the set of values ri — xi which have to be checked 

remain the same. Then max is attained at the end-point of this element r]^^'^^ because ri — xi is a 

n 

monotone increasing function of ri. 
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The last step in the proof is the final choice of ri. As was mentioned above ri must be an 
end-point of some element of r/^^^^ and takes finitely many values. Therefore ri should be chosen 

so that X2/^/N takes the largest possible value. Take the last point r[ = T^'^'i^^O on the orbit of 
of the length Qs^. Assume for definiteness that r[ lies to the left from 0. Consider m2 elements 
of 77^**^^ which start from r[ and go left. Then ri must be one of the end-points of these elements. 
Indeed, if ri lies more to the left from then the values xi take finitely many values and X2 will be 
significantly smaller. Therefore it cannot give maximum over r of our basic linear form. 

Thus we take m2 elements of ry^*^), consider their end-points. Each end-point is a possible value 
of r. Taking finitely many xi (see Lemma 1 and Lemma 2) we find minimum of our basic linear 
form. After that we find r for which this minimum takes maximal value. In this way we get the 
solution of our max-min problem. It is clear that this solution is a function of and elements 

hj, si < j < si + mi of the continued fraction of p near si. Since and hj, si < j < si + mi have 

limiting distribution as — > 00 the number /iv(a) = ^3^^ ^Iso has a limiting distribution. 

It remains to extend our proof to the case when the pairs from ai, 02, 03 have non-trivial common 
divisors, say ki is gcd of 01,03 and /c2 is gcd of a2,a3. It is easy to show that ki,k2 have a joint 
limiting probability distribution in the whole ensemble Q^. Fixing ki,k2 we can write oi = kia'^, 
02 = ^202, 03 = kik2a'^ where a'i,a'^ are coprime, 03,03 are coprime and fci,/c2 are coprime. This 
implies that (a'i)~^( mod Og) exists and we can multiply both sides of (3) by (o'^)"^. This will give 

kixi + A;202 • {a'i)~^ ■ X2 = ri(mod 03) (9) 

where ri = r • (a'^)~^ (mod 03). Denote b = 02(0'^)"^. 
Then from (9) we have the linear form 

kixi + k2 bx2 = Ti ( mod 03) (10) 

which we can treat in the same way as before. 



3 Statistical properties of continued fractions 

Statistical properties of elements of continued fractions usually are identical for real numbers 
and for rationales with bounded denominators (see [8]-[9]). 

Let M be the set of integer matrices S = {qq') with determinant detS* = ±1 such that 
1 ^ Q ^ Q', ^ P ^ Q, 1 ^ P' ^ Q'. For real a G (0, 1) the fractions P/Q and P'/Q' with 
S = {q Q') will be consecutive convergents to a (distinct from a) if and only if 
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(see ^ lemma 1]). Moreover if a = [0; hi, /i2, . . .] then for some s ^ 1 

^ = [0;/ii,...,/i,_i], ^ = [0;hu...,hs], (11) 

— = 0; /is, . . . , /ii , — - — — = 0; hs+i, /is+2, • • • • 
Q' -Qa + F 

It means that the distribution of partial quotients hs-k, ■ ■ ■ , /is+fc depends on Gauss-Kuz'min statis- 
tics of fractions Q/Q' and {Q'a - P')/{-Qa + P). 

For real a, xi, X2, yi, 2/2 £ (0,1) denote by Nx^,x2,y\,y2{(^i^) the number of solutions of the 
following system of inequalities 

0< 5-^0)^x1, Q^X2(5', Q^yiR, R^y2Q\ (12) 

with variables P, P', Q, Q' such that S = (g g/ ) G M. Let 

N{R) = Nr,^^x2,yi,y2{R) = / ^xi,x2,yi,y2{'^^R)da 

Jo 

and 

G(2;i, 2:2,^1,^2) = < 



^ f log(l + X1X2) log ^ - Li2(-xiX2) j , if a;2 ^ yiy2; 
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L c(2) 



Li2(-xiyiy2), if X2 > ym, 



where Li2(-) is the dilogarithm 

log(l - t) 



dt. 



The next statement implies Theorem [2j 
Proposition 1. For 2 



NiR) = G{xi,X2,yuy2) + 



xi log R 
R 



Proof. For every number a = [0;ai,a2, . . .] find a unique matrix S £ M with elements P, F' , Q, 
Q' defined by pT]) with the additional restriction Q ^ R < Q' . The inequalities < S^^{a) ^ xi 
define the interval Ixi{S) C (0, 1) of the length 



14,(5)1 



P' + xiF P' 



Q' + xiQ Q' 



Xl 



Q'{Q' + xiQ) 



Hence 



N{R) 



^ [Q^X2Q',Q^yiR,R^y2Q']Qr^Q?^y 
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where [A] is 1 if the statement A is true, and it is otherwise. Second row (Q, Q') can be comple- 
mented to the matrix from M in two ways. That is why 

Q'>R/y2 (Q,Q')=i ^ ^ 

In the first case X2 ^ yi2/2 and the Mobius inversion formula gives 

d^i? Q';}yiR/{x2d)Q^yiR/d^ ' 



log(l +2:1x2) log + / 

y X2 Jl/(xvx 



log(l + -^-l+0^"^'°^'' 



C(2) \^ X2 Jl/(xiX2) V t J \ R 

log(l + X1X2) log Li2(-xiX2j ] + ' 



C(2) V ^ " X2 -V . , ^ ^ 

The second case X2 > yiy2 can be treated in the same way. □ 

Let ^ 

^{R) ~ ^xi,x2,yi,y2{R) — ^xi,x2,yi,y2 [ti R 



(a,6) = l 



Theorem [3] will be proved in the following form. 
Proposition 2. For R ^ 2 



-L(i2) = G[xi,X2,yi,y2) +0 ' 



^K^-'J ^ 17 -^^7 nil a^j \ R 

Proof. Let a = a/& be a given number and S = (^q q,) £ M. he a, solution of the system (fT2]) . 
Denote by m and n the integers such that mP + nP' = a, mQ + nQ' = b. Then the system (I12p can 
be written as follows 

mP + nP' = a, mQ + nQ' = b, 
< m/n ^xi, < Q/Q' ^ X2, Q ^ yi-R, ^ y2(5'- 

Summing up solutions of this system over a and b we get that the sum L(R) equals to the number 
of solutions of the following system 

mQ + nQ' ^ i?^ < m/n ^ xi, < Q/Q' ^ X2, Q/yi ^ i? < y2Q', 
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where {qq') gA^, O^m^n, {m, n) = 1. For given Q and Q' values of P and P' can be founded 
in two ways. Number of solutions of the last system is equal to the area of the corresponding region 
with the factor 1/C(2) (see [HI Ch. II, problems 21-22]) 

+ 0' 



2C(2) Q'{Q' + xiQ) \ Q' 
It leads to the sum similar to (1131): 



Therefore 

^(^) = ^^(^) + 0(^1^' log' i?), 

and Proposition [2] follows from Proposition [TJ □ 

In order to prove theorem 4 we have to use Kloosterman sums 

Q 

Kq{m,n)= 2^ 1) e -j , 

x,j/=l 



where 5g(a) is characteristic function of divisibility by q: 

Sg{a) = [a = (mod p)] - 



1, if a = (mod 
0, if a ^ (mod 



Using Estermann bound (see [2]) 

\Kg{m,n)\ ^ ao{q) ■ {m,n,qf/'^ ■ q^l'^ . 
it is easy to prove the following statement (see [9] for details). 

Lemma 3. Let q ^ 1 be an integer, Qi, Q2, Pi, P2 be real numbers and ^ Pi,P2 ^ Q- Then 
the sum 

^g{Qi,Q2;Pi,P2)= '^9(^^-1) 



Qi<us;Qi+Pi 



satisfies the asymptotic formula 



$g(Ql, Q2; Pi, P2) = ^ • P1P2 + O ii^iq)) , 

where 

^l;{q)=ao{q) log\q + l)q'/^ . 
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It implies the following general result (see [8]). 

Lemma 4. Let q ^ 1 be an integer and let a{u, v) be a function defined on the set of integral 
points {u,v) such that 1 ^ u,v ^ q. Assume that this function satisfies the inequalities 

a{u,v) ^ 0, Aifia{u,v) ^ 0, Ao,ia(ii, t;) ^ 0, Ai_ia(ii, v) ^ (14) 

at all points at which these conditions have the well-defined meaning. Then the sum 

q 

W= ^ dq{uv — l)a{u,v) 

u,v=l 

satisfies the asymptotics 

^ u,v=l 

where ip{q) is the function from lemma\^ and A = a(l, 1) is the maximum of the function a{u, v). 
Let 



.b 



(a,i)) = l 



The next statement implies Theorem [H 
Proposition 3. For 2 



N,{R) =z-G{xi,X2,yi,y2)+0 



x\ log^ R 
RV^ 



^^(i?) =z ■ G{xi,x2,yi,y2) +0 



Proof. Let 

For a given z there is at most one matrix S = (^q q,) G M. such that Q ^ R < Q' and z G Ix^ {S). 
Hence 

iV.(i?)= 2: [Q ^ X2Q', Q ^ y.R, R ^ y2Q'] + ^ (|^) • 
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If Q' is fixed then P' and Q satisfy the congruence P' Q = ±1 (mod Q'). Therefore 
Q' 

Xl ( X\ 



E E '5QK^"Q±l)[Q^min{x2Q',yii?},P'^zQ']-— ^^ + 0(-^ 
Using Lemma m we obtain 

Q' 



'PiQ')sr\n^ ■ i n' mi , ^ log^ 



Applying the formula 

" ^ ^ ^^^^ 

we get the same sum as in the proof of Proposition [TJ 

As in Proposition [2] the sum Lz{R) equals to the number of solutions of the system 

mQ + nQ'i^R^, mP + nP' s^, z{mQ + nQ'), 
< m/n ^ Xl, < Q/Q' < X2, Q/yi ^ R< y2Q' , 

where {qq') S Ai, ^ m ^ n, {m,n) = 1. Again, there is at most one matrix 5* = {qq') 
such that Q ^ R < Q' and z G IxAS). Also for Q' ^ R 



This estimate implies that 



Lz{R) =^ E [^/y2 ^Q'^R^Q^ mm{y,R, X2Q'}] ^, ^^^^ ^^^^ + 0(a;ii?=^ log^ R) 



E E [Q ^ minj^ifi, ^ zQ'] t!-' + 0{x,R' log' R)- 



Using Lemma m one more time we obtain 

Q' 

^.(^) -|) IP [« < p' < + o « 
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Applying formula (jlSp we get the same sum as in as in the proof of Proposition [TJ □ 
Remark 1. In the simplest case X2 = yi = ?/2 = 1 w;e have cumulative distribution function 

F{x) = Fix, 1, 1, 1) = -^Li2(-x), 

which is not equal to the Gaussian function log2(l + x). As x ^ function F{x) (with error terms 
in Propositions \^ and \^ decreases as a linear function F{x) ~ 2x/C,{2). This fact shows that the 
expectation of the partial quotient (defined by inequalities qg-i ^ R < qg) is equal to infinity. 

4 Concluding remarks 

The calculations done by one of the authors (A. Ustinov) shows that the density of the limiting 
distribution of -^("i'"^'"^ j^as the following simple form: 



Pit) 



0, ifte[0,V3]; 
l|(tV3arccos*±^^ + |^/t^logg5|), if t G [2, +00) 



This result will be published elsewhere. 



References 

[1] BOURGAIN J., Sinai Ya. G. Limiting behavior of large Frobenius numbers. — Uspekhi Mat. 
Nauk 62 (2007), no. 4(376), 77-90; translation in the Russian Math. Surveys 62 (2007), no. 4, 
713-725. 

[2] ESTERMANN T. On Kloosterman's sum. — Mathematika, 8 (1961), 83-86. 

[3] Selmer E.S., Beyer O. On the linear diophantine problem of Frobenius in three variables. — 
J. Reine Angewandte Math., 301 (1978), 161-170. 

[4] Sinai Ya. G. Topics in Ergodic Theory, Princeton University Press, Princeton, NJ, (1994), 218. 

[5] Sinai Ya. G., Ulcigrai C. Renewal-type limit theorem for Gauss map and continued frac- 
tions. — Ergodic Theory & Dynam. Sys., 28 (2008), 643-655. 

[6] Sinai Ya. G., Ulcigrai C. A limit theorem for Birkhoff sums of non-integrable functions over 
rotations. — Probabilistic and Geometric Structures in Dynamics ( Contemporary Mathematics, 
Eds. K. Burns, D. Dolgopyat and Ya. Pesin. AMS, Providence, RI, 2008. 



12 



[7] Sylvester J.J. Problem 7382. — Educational Times 37 (1884), 26; reprinted in: Mathematical 
questions with their solution, Educational Times (with additional papers and solutions) 41 
(1884), 21. 

[8] Ustinov A. V. On the statistical properties of finite continued fractions — Zap. Nauchn. Sem. 
S.-Peterhurg. Otdel. Mat. Inst. Steklov. (POMI) 322 (2005), Trudy po Teorii Chisel, 186-211; 
translation in J. Math. Sci. (N. Y.) 137 (2006), no. 2, 4722-4738. 

[9] Ustinov A. V. On the number of solutions of the congruence xy = I (mod q) under twice 
differentiabie curve — Algebra and Analysis 20: 5 (2008), 186-216. 

[10] Ustinov A. V. On the distribution of Probenius numbers with three arguments I. — Mat. Sb. 
(to appear). 

[11] Ustinov A. V. On the statistical properties of elements of continued fractions — Doklady of 
the Russian Academy of Sciences, (accepted for publication). 

[12] Vinogradov I. M. Elements of number theory. — Moscow: "Nauka", 1972. 



October 29, 2008 :gpp 



13 



